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ABSTRACT 

We investigate the thermal stability of nuclear burning on the accreting X-ray pul- 
sar GRO J1744-28. The neutron star's dipolar magnetic field is 3 x 10 11 G if persis- 
tent spin- up implies that the magnetospheric radius is less than the co-rotation radius. 
Bildsten earlier noted that magnetic fields this weak might not quench the vigorous 
convection sometimes associated with thermonuclear instabilities, so that a convective 
burning front can propagate around the star in a few seconds and rapidly release the 
accumulated nuclear energy. After inferring the properties of the neutron star, we 
study the thermal stability of hydrogen/helium burning and show that thermonuclear 
instabilities are unlikely causes of the hourly bursts seen at very high accretion rates. 
We then discuss how the stability of the thermonuclear burning depends on both the 
global accretion rate and the neutron star's magnetic field strength. We emphasize that 
the appearance of the instability (i.e., whether it looks like a Type I X-ray burst or a 
flare lasting a few minutes) will yield crucial information on the neutron star's surface 
magnetic field and the role of magnetic fields in convection. 

We suggest that a thermal instability in the accretion disk is the origin of the long 
(~ 300 days) outburst and that the recurrence time of these outbursts is > 50 years. We 
also discuss the nature of the binary and point out that a velocity measurement of the 
stellar companion (most likely a Roche-lobe filling giant with mx 17) will constrain 
the neutron star mass. 
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1. Introduction 

The transient accretion-powered pulsar (v = 2.14 Hz) GRO J1744-28 was discovered by the 
Burst and Transient Source Experiment {BATSE) aboard the Compton Gamma Ray Observatory 
(CGRO) during a day of rapid bursting (about twenty bursts per hour) on 1995 December 2 
(Kouveliotou et al. 1996a; Finger et al. 1996a). Finger et al. (1996a) measured the 11.8 day orbit 
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and found that GRO J1744-28 was spinning up at a rate u = (3.5-12.2) x lCT 12 s~ 2 between 
1995 December 15 and 1996 January 23. The pulsar subsequently settled into a regime of hourly 
bursting with burst durations of 2-7 seconds, as seen by the PCA instrument on the Rossi X-Ray 
Timing Explorer (RXTE) (Swank 1996; Giles et al. 1996) and the OSSE instrument aboard CGRO 
(Strickman et al. 1996). A qualitatively similar burst has been seen from another accreting pulsar 
(see our discussion of SMC X-l in §2), but recurrent bursts of this nature have never been observed. 

The bursting behavior has two obvious energy sources: accretion and thermonuclear burning. 
Lewin et al. (1996) presented the case for accretion-powered bursts via analogy to the Rapid Burster. 
The rapid bursts (one every few minutes) are responsible for up to 50% of the total time averaged 
luminosity above 20 keV and thus cannot have a thermonuclear origin (Kouveliotou et al. 1996a). 
The hourly bursts (over 3000 observed by BATSE [Kouveliotou et al. 1996c]) are responsible for 
a much smaller fraction of the total luminosity, at least above 20 keV. The average burst fluence 
was 7 x 10~ 7 ergcm~ 2 (20-50 keV) on 1996 January 15, when there were 40 per day and the 
persistent luminosity was 2.5 Crab in the 20-100 keV band (1996 January 16 [Fishman et al. 1996]) 
and 4.4 ± 0.3 Crab in the 8-20 keV band (1996 January 14-1996 January 15 [Sazonov & Sunyaev 
1996]). Above 20 keV, this gives a time-averaged burst luminosity 60-100 times smaller than the 
steady-state accretion luminosity. The burst energetics below 20 keV (where most of the energy 
is emitted) is still uncertain due to dead-time corrections in the PCA. Preliminary indications are 
that these corrections are large enough so that the bursts might be responsible for more than 5% 
of the < 20 keV emission (Jahoda 1996, private communication). 

The burst energies were marginally consistent with nuclear energy release and motivated our 
theoretical study of thermonuclear burning on this unusual X-ray pulsar. Additional motivations 
were the matchings of the characteristic decay time (2-10 seconds [Strickman et al. 1996]) with the 
cooling time at the nuclear burning depth and the mean recurrence time (about 30 minutes for the 
260 bursts seen by OSSE from 1996 January 16 to 1996 January 30 [Strickman et al. 1996]) with 
the time to accumulate enough fuel for an instability. The primary observational evidence against 
all of these bursts having a thermonuclear origin is (1) the independence of the recurrence time 
from the accretion rate (at least when it is brighter than 200-400 mCrab [Giles et al. 1996]), (2) 
the existence of "foreshocks" before the bursts (Giles et al. 1996), (3) the lack of spectral evolution 
during the burst, and (4) the global burst energetics from the PCA instrument (if the dead-time 
corrections are understood). 

We consider here the possibility that this X-ray pulsar has an unusually low field and therefore 
might exhibit different thermonuclear burning behavior than conventional X-ray pulsars. It is 
unlikely that the thermonuclear burning was unstable at the peak of the outburst. Hence we 
concentrate on lower accretion rates, for which the burning will most likely be thermally unstable 
and will manifest itself as Type I X-ray bursts or flares of a few minutes duration. 

Since the presence and character of a thermonuclear instability depends on the neutron star's 
magnetic field (B) and global accretion rate (M), we begin in §2 by summarizing the indirect 
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inferences about these quantities. The spin behavior points to an especially weak (<C 10 12 G) 
dipolar magnetic field component and constrains the global accretion rate as well. We also compare 
GRO J1744-28's bursting behavior at the peak of the outburst to that of SMC X-l, for which 
similar arguments point to a low dipole field. In §3, we constrain the binary properties, infer a 
time-averaged mass transfer rate, and speculate on the origin of the long-term transient behavior. 
The thermonuclear stability of the accreted hydrogen and helium is discussed in §4 for a range of 
M's. Section 5 is a discussion of the magnetic field's role in the stability and character of nuclear 
burning in accreting X-ray pulsars. In particular, we discuss in detail how the star will behave as 
M decreases. We conclude, in §6, by describing what a successful identification of a thermonuclear 
instability implies about the neutron star's properties. 

2. Properties of the Accreting Neutron Star and Comparison to SMC X-l 

In addition to the bursting behavior, this pulsar is unusual because of its unusually high 
spin frequency and its steady spin- up over a large range of accretion rates. In the context of 
magnetic accretion, these facts imply a dipole field lower than most other accreting pulsars. The 
magnetosphere is located at r m = ^r^ (Ghosh & Lamb 1979), where is a characteristic length 
found by equating magnetic and fluid stresses, and £ is a model-dependent dimensionless number. 
Estimates of £ range from ~ 0.52 (Ghosh & Lamb 1979) to ~ 1 (Arons 1993; Ostriker & Shu 1995; 
Wang 1995). A measurement of £ obtained from the observed quasi-periodic oscillations in the 
accreting pulsar A0535+26 (Finger, Wilson, & Harmon 1996b) gave £ ~ 1. Since the neutron star 
is spinning up, the magnetospheric radius, r m , must be less than the co-rotation radius, r co , which 
implies an upper limit on the surface strength of the dipolar componenet of the magnetic field 



in agreement with previous estimates (Finger et al. 1996a; Sturner & Dermer 1996). The magne- 



This upper bound is less than the typical X-ray pulsar field strength, which, as we discuss in §5, 
changes the nature of the convection; GRO J1744-28 is the first high accretion rate X-ray pulsar 
beneath this limit. 

The X-ray spectrum above 20 keV falls very steeply (Kouveliotou et al. 1996a; Strickman et 
al. 1996) and is consistent with being above the characteristic cut-off energy found by RXTE/PCA 
(pa 15-20 keV [Swank 1996; Giles et al. 1996]). Daumerie et al. (1996) argue that this spectrum 
and the increase in pulse fraction with energy imply a surface field ~ 10 12 G. This argument is at 
odds with our upper limit on the dipolar component (for £ = 1) and might imply that higher order 
magnetic moments are present. 

We use the measured spin-up rate i> to constrain the accretion rate onto the neutron star. 
The maximum specific angular momentum of the accreted matter is Z max = {GM x r co ) x l 2 , where 




(1) 



tospheric radius is presently unknown; continual spin-up at lower M's will reduce this upper limit. 
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r co = 1.0 x 10 8 cm ss lOOi? is the co-rotation radius (where the Kepler frequency equals the neutron 
star's spin frequency) for a M x = 1.4M neutron star. Because the observed torque lixlv must be 
less than M/ max (Ghosh & Lamb 1979; Chakrabarty et al. 1993), we set a lower bound to M, 

M>8x 1 0-»M oyr -( IF ^)( I A ;; ) 2 , (2) 

where / = 0AM x R 2 is the neutron star's moment of inertia (a good approximation for our choice 
of mass and radius [Ravenhall & Pethick 1994]). The spherically averaged local accretion rate, 

™ sph = 4^ > 4 X g Cm " 2 ( lO-n s -0 ' (3) 

is then independent of the neutron star's radius. Using a rough bolometric flux for 1996 January 
15 of 10~ 7 ergcm~ 2 s _1 , we infer a minimum distance d = 3kpc. As we note later, the extinction 
and position in the galaxy most likely place the object at least 8 kpc away. 

The only other accreting pulsar which has always been spinning up (z> = 2.4 x 10~ n s~ 2 ) and 
for which our earlier arguments yield a comparable field strength is SMC X-l (y = 1.4 Hz). For a 
typical accretion rate M ss 4 x 10~ 8 M Q yr -1 (Levine et al. 1993), we infer B < 10 12 G. Angelini, 
Stella, and White (1991) discovered an X-ray burst from SMC X-l during an EXOSAT observation 
on 18 October 1984 when the 1-16 keV luminosity was L ~ 4 x 10 38 ergs -1 . This burst is very 
similar to the large hourly bursts seen from GRO J1744-28. The burst rose by a factor of three 
within one second, lasted for about 80 seconds, and was followed by a 35% decline in the persistent 
flux. There was no evidence for spectral changes during or after the burst, and the pulse fraction 
and phase remained constant. The recurrence time must be long, as only one burst was seen in 
~ 20 hours of observation. Angelini et al. (1991) noted the similarities to the Rapid Burster and 
argued for accretion as the energy source for these events. They also noted that the variability of 
the source has a strong underlying quasi-period of a few minutes, which was present in all EXOSAT 
observations except for the five hours following the burst. About 10% of the total luminosity of 
the source is in these variations. 



3. Properties of the Binary 

Following the ROSAT positioning of this object (Kouveliotou et al. 1996b), Augusteijn et 
al. (1996) identified the variable infra-red counterpart, which was present in an earlier image (8 
February 1996) of Blanco, Lidman, and Glazebrook (1996) at mx = 15.7 ±0.3 and was undetected 
and at least a magnitude fainter on 28 March 1996. This light is most likely X-rays reprocessed 
by either the accretion disk or companion (Augusteijn et al. 1996). The most obvious Roche- lobe 
filling object is a first-ascent red giant branch star with a degenerate helium core of mass Mh c 
and an overlying hydrogen envelope (Finger et al. 1996a; Sturner & Dermer 1996). Hydrogen shell 
burning via the CNO cycle supplies a luminosity strongly dependent on only the helium core mass, 
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which allows us to estimate the stellar luminosity and expected IR magnitude. In the absence of 
mass or angular momentum loss, this binary evolves due to the expansion of the red giant as the 
helium core mass grows. 

3.1. Constraints on the Optical Companion and Neutron Star Mass 

We solve for the companion mass, M c , by using the orbital parameters P or b = 11.83 d and 
a x sini = 2.631t-sec (Finger et al. 1996a) and the core-mass/luminosity relations of Webbink, 
Rappaport, & Savonije (1983). We presume that the giant fills the Roche lobe estimated by 
Eggleton (1983) and fix M x = 1.4M . For metallicities Z = 0.02(0.0001), the first allowed solution 
(corresponding to an unrealistic hydrogen envelope mass of zero) has M c /M = 0.216 (0.232) and 
inclinations less than 18-19 degrees (cosi > 0.95). The solution in the middle of the allowed range 
(i.e., cosi = 0.975) has M c = 0.334 M , M He = 0.22 (O.24)M , L « 12L , an orbital separation of 
26i? , a stellar radius of 6.9i? , and a projected companion velocity K ~ 20kms~ 1 . The helium 
core mass is slightly larger for the metal-poor case to compensate for fewer catalysts. The strong 
dependence of the giant radius on Mh c makes the inferred core mass (shown by the dotted line 
in Figure 1) nearly independent of the inclination angle and always close to 0.22 (O.24)M . The 
dashed line in Figure 1 shows M c as a function of cos i for Z = 0.02 and M x = 1.4M . The resulting 
mass transfer rate (Webbink et al. 1983) is given by the solid line in Figure 1. For the "typical" 
solution presented above, the average mass transfer rate is M ss lO~ 9 M yr -1 , and the hydrogen 
envelope mass is O.1M . This implies a lifetime, nearly independent of the metallicity, of 10 8 yr. 
Most of the giant's envelope goes onto the neutron star (for conservative evolution). 

Unless we are looking at the system nearly pole-on, we must conclude that M c ss O.3M . If 
this low-mass star followed a normal evolutionary track, then it must have started mass transfer as 
a ~ Mq star. If all the departed mass accreted onto the compact object, then either the neutron 
star is more massive than its "birthweight" (1.2-1.7M [Timmes, Woosley, &: Weaver 1996]), or 
the neutron star was a white dwarf that underwent accretion-induced collapse. Considering more 
massive neutron stars does not alleviate the companion's substantial mass loss but does increase 
the inclination angle (e.g., cosi = 0.9 for M x = 2.OM , a doubling of the allowed phase space). 

Depending on the inclination angle, an infrared measurement of the orbit of this system might 
actually constrain the neutron star mass. The projected companion velocity is 



We find the minimum companion mass by requiring (1) that the star be luminous enough so that 
it fills the Roche lobe, and (2) that the overlying hydrogen envelope have a minimum mass of 
1O _2 M so as to live for at least 10 7 yr and to be fully convective. The minimum companion mass 
(M c = O.22M ) is then basically independent of M x , so that the maximum K for the Roche- lobe 
filling giant hypothesis is K mSLX = 2\~km.s~ 1 (M x /M & ). Hence, a measurement of K can potentially 




(4) 
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constrain the neutron star mass. For example, if M x = 2M & and we presume that cos i is uniformly 
distributed between 0.9 and 1.0, then there is a 50% chance of measuring a velocity larger than 
the maximum for the 1.4M Q neutron star case. This will be a tough job. We find that the stellar 
effective temperature is in the range 4000-4200 K, which, for the luminosity derived above, gives 
M K ~ -0.2 (Bessell & Brett 1988). If the extinction is Ay k, 25, as implied by the N H ~ (5- 
6) x 10 22 cm" 2 measurement of Dotani et al. (1996), then Ak ~ 3 and for a distance of 8 kpc, we 
expect mx ~ 17. 

3.2. The Long-Term Accretion History 

Consistent with the transient nature of this binary, the inferred long-term mass transfer rate is 
less than the current accretion rate. This possibly indicates that the present outburst (which so far 
has lasted for over 300 days) arises from a thermal instability in the accretion disk similar to that 
occurring in dwarf novae. Indeed, the implied long term mass transfer rate is much less than the 
amount needed for stable mass transfer (Shatter 1992). Our orbital calculations give the distance 
from the neutron star to where the accreting matter strikes the accretion disk as R r ~ 4i? Q (Lubow 
& Shu 1975). The thermal instability occurs when the column density of matter (£) accumulated 
at this radius exceeds ^ 10 3 gcm -2 (Cannizzo & Wheeler 1984), in which case the accumulated 
mass is ;> i? 2 S ~ 5 x 10~ 8 Mq. This requires over 50 years of mass transfer from the giant and 
supplies more than enough material to power the large outburst. 

The accretion will most likely be halted by a propeller effect at low accretion rates (Illarionov 
& Sunyaev 1975). However, even when not accreting, the neutron star has a luminosity from the 
cooling core (at temperature T c ) given by Gudmundsson, Pethick, and Epstein (1983) as 

This luminosity will probably not quench the thermal instability of the accretion disk as fully 
exposed matter at R r only has an effective temperature of T c g = 1800K(T C /10 8 K) a55 . It is the 
unusual combination of large orbital separation, low inferred mass transfer rates, and low neutron 
star luminosity that allows for such instabilities in this neutron star binary. For time-averaged 
accretion rates ~ 10 -9 MQyr -1 , the outburst intervals are longer than the history of X-ray and 
7-ray monitoring. The system's extinction (Ay > 25) also makes optical identification of a prior 
outburst unlikely. 

4. The Nuclear Burning of Accreted Matter 

The thermal stability and appearance of nuclear burning on steadily accreting neutron stars is 
well studied. For comparable metallicities and magnetic fields weak enough so as not to affect the 
opacities (< 10 13 G), the only residual difference between this neutron star and others accreting at 
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comparable local rates is a colder core. This might allow for unstable hydrogen/helium ignition at 
slightly higher instantaneous local accretion rates (m) than on a steadily accreting star. 

Constant accretion at these M's conductively heats the neutron star core, which attains an 
equilibrium temperature of T c = (2-4) x 10 8 K in about 10 3 -10 4 yr by balancing this heating with 
neutrino cooling (Ayasli & Joss 1982). The time-averaged M in GRO J1744-28 is comparable to 
many X-ray binaries, whereas the instantaneous M can clearly be much higher. The cycling of 
M leads to a colder equilibrium core than that of a neutron star steadily accreting at the same 
time-averaged M. This is because the thermal timescale of the deep ocean and crust (which is 
where substantial energy release occurs as matter is forced through the electron capture boundaries 
[Haensel & Zdunik 1990]) is ~ 10-100 years, much longer than the outburst duration. There 
is thus insufficient time to heat the crust to a temperature profile favorable for sending a large 
luminosity L « (M/m p )(l MeV) into the core. In addition, between outbursts, the cold outer 
envelope conducts away the accumulated heat in the deep ocean and crust. For temperatures 
below T c ~ 2 x 10 8 K the core cools radiatively between outbursts, much as a young neutron star 
still hot from birth. 



4.1. Settling and Nuclear Burning of the Accreting Matter 

As described in §3, the neutron star is accreting the hydrogen-rich envelope of an evolved 
giant. Because the giant has already lost an appreciable amount of mass, the matter presently 
being transferred was most likely processed in the companion's interior during its main sequence 
lifetime. This would increase the helium content. The /3-decay limitations and high temperatures 
on the neutron star fix the hydrogen burning (via the CNO cycle) rate at the /3-decay limited value 
of the "hot" CNO cycle, en = 5.8 x 10 15 Zcno erg g -1 s _1 where ^cno is the mass fraction of the 
CNO elements. For high accretion rates, this burning never consumes the accreted hydrogen before 
helium ignition, so that unstable helium burning occurs in a hydrogen-rich environment, which 
enhances the nuclear reaction chains and energy release (Lamb & Lamb 1978; Taam & Picklum 
1979; Fujimoto et al. 1981; Taam 1982). Because of the complicated rp process (Wallace & Woosley 
1981), the exact composition of the ashes is unknown except for a few limited cases (Van Wormer 
et al. 1994). 

Prior to helium ignition, the accreted material is settling onto the neutron star and is stably 
burning hydrogen. Since the thermal time is always less than the time to accrete to a given depth, 
^accr = y/ih, we find the temperature structure by solving the time- independent entropy and flux 
equations, 

m ds dF _ c daT 4 , , 

where the local accretion rate is written as m = m 4 10 4 g cm~ 2 s _1 , s is the specific entropy, F 
is the outward heat flux, and y is the column depth. The opacity in the upper atmosphere is 
given by 1/k = l/ft r ad + 1/ftcond) where the radiative opacity is the sum of electron scattering (we 
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usc Paczyhski's [1983] fitting formulae for the degeneracy and high temperature corrections) and 
free-free absorption. We use the conductivity from Yakovlev & Urpin (1980). 

These settling solutions are valid until the helium burns fast enough to appreciably change 
either the temperature or the helium abundance. Fushiki and Lamb (1987) defined the boundary 
of stable helium burning in the y-T plane (the ignition curve) by setting 

dT ~ cff' U 

where e c = acT 4 /3ny 2 is a local representation of the conductive cooling and e^ a is the helium 
burning rate. These derivatives are taken at constant pressure, as the instability grows faster than 
the pressure changes. The heavy dot-dashed line in Figure 2 denotes this ignition curve for Y = 0.3. 
Regions to the right of this curve are thermally unstable to helium burning. We also define the 
depletion curve by equating the helium lifetime to the 3a reaction with i accr (Fushiki & Lamb 
1987); the heavy solid lines in Figure 2 denote this condition for rh^ = 3, 7.5, and 30. 

We derive the settling solutions by varying the flux exiting the atmosphere until the flux is 
zero at the depth where the ignition curve (or depletion curve, whichever is first) is met. We take 
zero flux at the bottom to mimic the cold core. These settling solutions all strike the ignition curve 
before the depletion curve, so that we find the column density 7/ ign accumulated on the star prior 
to unstable helium ignition. 

For Population I type metallicities, the temperature of the settling material is mostly set by 
the slight amount of hydrogen burning. Indeed, as is evident in the three light solid lines that 
display these settling solutions in Figure 2, y- l&i depends only mildly on rh for the Zqno = 0.01 
case. We obtain y ign /10 8 g cm~ 2 = 3.1, 2.95, and 2.65 for m± = 3, 7.5, and 30; these correspond 
to recurrence times of 2.9 hours, 1.1 hours, and 15 minutes. The flux exiting the atmosphere is 
F/(10 22 ergcm- 2 s _1 ) = 2.0, 2.42, and 4.93 for these m's. 

Lower CNO abundances are relevant if the companion is an older Population II giant and/or 
if spallation of the incident nuclei occurs in the accretion shock (Bildsten, Salpeter, & Wasserman 
1992). At lower Zcno> the settling solutions are more sensitive to the accretion rate via the 
gravitational compression terms. The two bottom thin solid lines in Figure 2 are settling solutions 
when Zcno = 10~ 4 for 777,4 = 7.5 and 30, giving recurrence times of 6.6 hours and 0.5 hours. 

4.2. When is the Burning Time Dependent? 

The crucial question to answer is, "At what accretion rate is the burning unstable?". The 
simplest way to address this is to construct the steady-state solution (represented by the dotted 
lines in Figure 2 for the m's given above) that burns the material as fast as it accretes. If this 
solution does not consume the fuel before reaching the instability curve, then it is unstable. For an 
accretion rate below m cr i t ~ 3 x 10 4 gem -2 s -1 , the burning is unstable (for Y = 0.3). This critical 
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accretion rate increases by a factor of two if the helium mass fraction is 0.5. These estimates also 
indicate that the burning is stable when m > w. C rit-[| Since m cv \t is smaller than our minimum 
estimated accretion rate at the outburst peak, it is unlikely that the hourly bursts during the brightest 
parts of the outburst have a thermonuclear origin. 

However, the burning becomes unstable when m < m cr ;t (either due to a reduction in the 
overall accretion rate or spreading away from the polar cap; see §5). When the fuel spreads over 
the whole star prior to ignition, the maximum luminosity (and hence flux) at which an instability 
occurs is L ~ 7 x 10 37 ergs -1 (for m cr it = 3 x 10 4 gcm~ 2 s _1 and R = 10km). This corresponds to 
a bolometric flux 

^unstable < KT 8 erg cm"' ( 7x J eigs ^ ) (^p) 2 . (8) 

We convert this to cs _1 in the PC A instrument aboard RXTE by using the conversion factor 
(Giles et al. 1996) of 4 x 10" 8 erg cm" 2 s" 1 (2-60 keV) for 10 4 cs _1 . For d = 8kpc, unstable 
burning signatures will most likely not appear until the PCA count rate is <^ 2500 cs" 1 . 

When the burning is unstable, the intersection of the settling solutions with the helium ignition 
curve defines 7/i gn . For typical metallicities, yi gn > 1.5 x 10 8 gem -2 , which at m C{ \t gives a recurrence 
time of just over an hour. Once the helium ignites, the temperature rises rapidly and starts a 
propagating combustion front. If the whole star is ignited (see next section) the maximum burst 
energy would be 47ri? 2 yig n (7MeV/m p ). This is the maximum, as time-dependent calculations of 
the hydrogen/helium burning flash often found incomplete burning (Taam et al. 1993) and the 
whole star need not ignite. The resulting fluence would be 

Maximum Burst Fluence = 1.5 x 10 -6 erg cm -2 ( ^ ( R . (9) 

& \1.5 x 10 8 gcm-V VlOkm/ \ d J KJ 

The peak luminosity depends on the combustion front's propagation speed through the fuel-rich 
regions, as we now discuss. 

5. The Role of the Magnetic Field 

Conventional Type I X-ray bursts are not seen from highly magnetized (B ;> 10 12 G) accreting 
X-ray pulsars. This was at first surprising because they accrete at rates comparable to X-ray 
burst sources that are not obviously magnetic. Joss and Li (1980) explained the lack of bursts 



1 There have been indications in the past that the rh C rit obtained from a time dependent simulation might actually 
be higher than our estimate. In particular, Ayasli and Joss (1982) argued that the steady state solution would not 
be reached until the time to reach the ignition depth became shorter than the local thermal time. This requires 
accretion rates 2-5 times higher than our estimate. However, in the absence of a full time-dependent calculation, we 
will stick to our present, and maybe overly conservative, estimate. 
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by stabilizing the nuclear burning with two mechanisms. For this pulsar, the relevant mechanism 
is the increased local accretion rate on the polar cap. The magnetic field funnels the accretion 
onto the polar cap and confines the accretion mound until the ignition pressure is reached. These 
constraints are typically satisfied for steadily accreting X-ray pulsars with M ^ lO -lo M0yr -1 , as 
the fractional area of the polar cap only needs to satisfy A cap /47ri? 2 <^ 0.01. This is well within the 
estimates obtained by either following the field lines from the magnetospheric radius to the star 
(Lamb, Pethick, & Pines 1973), 

M \ 2/7 f R \ 9 / 7 
lO -8 M yr -1 J v 10 km J ' (10) 

or allowing the matter to penetrate through the magnetopause via a Rayleigh- Taylor instability 
and attach to field lines at smaller radii (Arons & Lea 1976; Eisner <fe Lamb 1977). 

The absence of Type I X-ray bursts in all X-ray pulsars seems unlikely given the range in 
magnetic field strengths and accretion rates. Bildsten (1995) suggested that, even when the burning 
is thermally unstable, a strong magnetic field will inhibit the rapid lateral convective motion needed 
for the combustion front to ignite the whole star in a few seconds (Fryxell & Woosley 1982). The 
field strength required to halt the convective (~ 10 6 cms -1 ) propagation of burning fronts is not 
known. Convection is potentially stabilized when B 2 > 8ttP (Gough & Tayler 1966), which requires 
B ^ 7 x 10 11 G in the helium burning region, f] Most inferred dipolar field strengths for accreting 
X-ray pulsars easily satisfy this constraint. However, GRO J1744-28 does not and hence is an 
especially intriguing candidate for showing Type I X-ray bursts. 

If the lateral convective motion is inhibited (or if there is too little fuel for convection to 
occur [Bildsten 1993]) then the burning front propagates at the slower speed set by heat transport 
(electron conduction and/or radiative transport, depending on the depth) 

Vslow ~ 80 + 200 ( - m °~ Vq cms" 1 , (11) 

\4.5 x 10' gcm -i / 

where yjje is the local helium column density (in gem -2 ) and y q = 1.35 x 10 s gem -2 is the min- 
imum column density needed for a pure helium burning front to propagate (Bildsten 1995). This 
relation is for a pure helium atmosphere and is a reasonable lower limit to the mixed hydro- 
gen/helium burning case. For the typical yne ~ 2 x 10 8 gcm -2 where the instability sets in, 
v s \ ow ~ 400 cms" 1 , so that the burning front crosses a 10 5 cm polar cap in i C ross ~ 4min. In this 
case, the thermonuclear instability would appear as a flare of a few minutes duration (potentially 
time symmetric) with a luminosity set by the amount of accumulated fuel and the time to burn all 
of it, L flare sa 0.05L accr (t accr /t cross ) (Bildsten 1995). 



2 Even lower fields might slow the burning fronts, as the sub-sonic velocities (v c ~ 10 6 cms -1 ) implied by efficient 
convection can only push around fields of strength B 2 < Sivpv 2 , or B < 10 9 G at the helium ignition depth (Bildsten 
f995). Observations will most likely tell us the outcome for magnetic fields in the intriguing regime pv 2 <C B 2 /8tv <C 
P(10 9 G < B < 7 x 10 11 G). 
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5.1. The Spreading of Accreted Fuel 

The local accretion rate, m, onto the polar cap would greatly exceed m sp h if accretion onto 
GRO J1744-28 is via a polar cap of the size given by equation (^) (this is not clear given the 
small pulse fraction and nearly sinusoidal pulses seen by Finger et al. [1996a]). Because the fate of 
the accreted material depends on m, knowing the depth at which accreted matter flows laterally 
over the surface and reduces m to m sp h is crucial. The magnetic Reynolds number of the flow is 
Anaiv/c 2 , where £ is the length over which the magnetic field varies, a is the electrical conductivity, 
and v = m/p is the downward flow velocity. Using the pressure scale height, h = P/pg, as an 
estimate for £, we find that the magnetic Reynolds number is ~ 100 in the helium burning region 
(P ~ 10 22 erg cm -3 , T ~ 5 x 10 8 K). As a result, the magnetic field is frozen into the matter at the 
ignition depth. 

For a small polar cap, we find the pressure where the magnetic field can no longer hold up 
the accretion mound by balancing the transverse pressure gradient (given by the lateral extent of 
the polar cap, -^cap ^ \/ Acap ^ P) with a distorted field. Following Hameury et al. (1983), we 
presume spreading occurs at the depth where the field is distorted by a large angle. To illustrate 
this problem, consider an azimuthally symmetrical and poloidal magnetic field B = (B w ,0, B z ) 
in cylindrical coordinates (zu,(j),z). The characteristic length in the z-direction is the pressure 
scale height h = P/pg and the characteristic length in the ro-direction is £ C ap- The accretion flow 
distorts the field from B = (0, 0, B) to the perturbed configuration B = (P ro , 0, B — 5B Z ). Equating 
estimates of J obtained from curl_B = AirJ/c and the ro-component of —VP + pg + JxB/c = 0, 
and using div_B = to obtain a relation between 5B Z and B m , we have 

P m AirhP 



B 4apP 2 ' 



(12) 



For a fully ionized H/He mixture with an ideal gas equation of state, the pressure where P ro ~ B 
is 

w ^cap \ 1/2 [ B \ 2 /2xl0 8 K\ 

~ 10 ergcm (io^) [whO [—t-)- <13> 

This pressure defines a boundary {y s in column depth) where the matter starts to spread laterally. 
We will presume that above y s the local accretion rate is M/A cap and that below y s , the local 
accretion rate is the spherical value m sp h- 



5.2. Global Behavior of the Nuclear Burning for GRO J1744-28 

Having shown that the nuclear burning is unstable when fn < m cvlt ~ 3 x 10 4 gcm _2 s _1 , we 
now estimate rh at the ignition depth in terms of the dipole magnetic moment p and global accretion 
rate M (we are assuming for simplicity that no higher-order multipole moments are present). We 
show in Figure || the cases £ = 1.0 and £ = 0.5. Within this parameter space, the first requirement 
is that the magnetospheric radius be less than the co-rotation radius (eq. [Q] with R set to 10 km) 
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and is indicated by the unhatched area. The second relation, indicated by the dot-dashed line in 
Figure ||, is y s = y\ zn . Below this line, the accretion flow spreads before ignition, thereby lowering 
the local accretion rate to m S ph- Depending on the dipolar field strength there are two cases. 

• For B <^ (2-4) x 10 10 G, the accreted matter spreads before igniting regardless of the polar 
cap area. Thus, the relevant local accretion rate (indicated by the vertical dashed line) is 
m sp h and the burning is unstable when the flux constraint (eq. ||) is reached. This case is 
denoted by region I in Figure 3. 

• For B ^ (2-4) x 10 10 G, the matter ignites prior to spreading, so that the nature of the 
thermonuclear burning depends on the polar cap size. This case is denoted by region II in 



Figure 3. The most constraining polar cap area is given by equation 10 and is indicated by the 
slanted dotted line in the right-hand plot. For GRO J1744-28, this scenario is only relevant 
if £ ~ 0.5, as otherwise the propeller effect will halt accretion before m < m cr ;t on the polar 
cap. For £ = 1.0, nuclear burning in region II is unstable only if the polar cap area is larger 
than the estimate of equation (|To|). 

Our ignorance of the polar cap size prohibits us from saying at what accretion rate the B ;> 
(2-4) x 10 10 G case becomes unstable (the triangular regions denoted II). If we choose a fixed polar 
cap area of 10% of the stellar area, then the M needed for an instability decreases by a factor of 
ten, as does the maximum flux required to see an instability. The polar cap of Arons and Lea 
(1980) is always larger than 10% of the stellar area and implies that the burning becomes unstable 
when L < 3 x 10 37 ergs~ 1 . The appearance of unstable burning will constrain the polar cap size, 
as a necessary condition for stable burning is A ca _ p < M/m cr i t . 



6. Summary and Observational Outlook 

The continuous monitoring of GRO J1744-28 by the RXTE provides an important opportunity 
to learn about both accretion and thermonuclear instabilities on a weakly magnetized neutron star. 
We have shown that the bursts observed during the peak of the outburst are most likely not of 
thermonuclear origin, as even the minimum local accretion rate on the neutron star is too high for 
unstable burning of hydrogen-rich material. This statement is no longer true as the accretion rate, 
M, decreases. 

The full understanding of the thermonuclear burning depends on many properties of both the 
neutron star and the binary. In §2, we used the observed torque to estimate the minimum accretion 
rate and dipole field, which led to the comparison with SMC X-l, another bursting X-ray pulsar. 
We speculated in §3.2 that the present outburst might be the result of a thermal instability in the 
disk, comparable to what occurs in dwarf novae. We also discussed the optical companion and 
pointed out that a careful velocity measurement of the companion will constrain the neutron star 
mass. 
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We have shown that, if B <^ (2-4) x 10 10 G, then the nuclear burning becomes unstable when 
the intrinsic source luminosity is L <^ 7 x 10 37 ergs~ 1 . For higher magnetic fields, the matter stays 
confined at the polar cap before ignition, in which case the stability of the nuclear burning depends 
strongly on the polar cap size. As the outburst fades and M decreases, the burning might become 
unstable, especially if the polar cap is larger than the conventional estimate (equation ]lQ|), for 
which the thermonuclear burning is always stable when £ = 1. Hence, for larger fields, the global 
accretion rate when the first signatures of thermonuclear instability appear will constrain the polar 
cap size. 

As discussed in §5, the character of the instability strongly depends on the burning front's 
propagation speed. If the instability is too weak to convect (or if the field is strong enough, 
B ^ 7 x 10 11 G, to inhibit the lateral convective heat transport) a possible burning signature is 
flares with durations of a few minutes to an hour. The duration of the flare is set by the slow 
burning speed (§5) and the size of the fuel-rich region. These flares rise on a timescale comparable 
to their duration and are not necessarily asymmetrical in time, as Type I X-ray bursts are. For a 
polar cap radius of 10 5 cm (10% of the stellar radius) and a typical ignition depth of 2 x 10 8 gem -2 , 
the burning front crosses the polar cap in four minutes. If convection occurs, then we expect to see 
Type I X-ray bursts as in other Low Mass X-Ray Binaries. 

When bursts or flares first appear, they will have recurrence times of order an hour if our 
m cr it determination is accurate. It is possible that the recurrence times will be longer than this 
(especially for low metallicity), in which case detection will be more difficult. The ability of the 
PCA to position bursts to within 0.2 degrees has already eliminated GRO J1744-28 as the source 
of two Type I X-ray bursts seen in the field (Corbet & Jahoda 1996; Jahoda et al. 1996) and should 
eventually provide an unambiguous localization of a Type I burst from GRO J1744-28. 

We thank Jon Arons, Mark Finger, Keith Jahoda, Chris McKee, Ed Morgan, Tom Prince 
and Bob Rutledge for many discussions about the nature of this source. We especially thank Bob 
Rutledge (MIT) for creating and maintaining a Web site about this X-ray binary. Our work was 
supported by NASA via grants NAG 5-2819 and NAGW-4517 and by the California Space Institute 
(CS-24-95). L. B. was also supported by the Alfred P. Sloan Foundation. 
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Fig. 1. — The companion mass and mass transfer rates as a function of inclination angle for 
M x = 1.4M . For the orbital parameters of Finger et al. (1996), we show the red giant companion 
mass (dashed line) and inferred mass transfer rate (solid line) as a function of cosi. There are no 
Roche-lobe filling solutions for i > 18 degrees if M x = 1.4M Q . 
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Fig. 2. — Helium ignition for hydrogen-rich accretion (Y = 0.3) at accretion rates appropriate for 
this source for a 1AM@, R = 10km neutron star. The heavy solid lines show the helium depletion 
curves for 777,4 = 3, 7.5, and 30. The heavy dot-dashed curve is the helium ignition curve. The top 
three thin solid curves are settling solutions for Zcno = ^-01 an d rhi = 3, 7.5, and 30, while the 
bottom two are for Zcno = 10" 4 and 7774 = 7.5 and 30. The three dashed lines are the envelope 
structure if the burning were in steady-state for these three m's. 
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Fig. 3. — A graphical summary of the burning behavior parameterized by the dipole magnetic 
moment fi and the global accretion rate M. The left plot is for £ = 1.0; the right, for £ = 0.5 
(cf. equations and [[Hj). We assume for simplicity that no higher-order multipole moments 
are present. In the shaded region the propeller effect stops accretion (i.e., r m > r co ). For the 
region beneath the dot-dashed line labeled y s = yi gn , spreading occurs above the ignition column 
depth «2x 10 8 gem -2 . The vertical dashed line indicates the critical accretion rate, m cr i t = 
3xl0 4 gem -2 s _1 , for spherically symmetric accretion. In region I, the burning is unstable regardless 
of polar cap size. For region II, the accreted material is magnetically confined at the ignition 
column density; our ignorance of the polar cap size prohibits us from indicating where the burning 
is unstable. In the right-hand plot we show with a dotted line where unstable burning begins if 
the accreted material is confined to a polar cap of area given by equation (10). Nuclear burning is 
unstable in the region to the left of this line. 



